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ABSTRACT 

Skin  damage  thresholds  were  measured  and  compared  with  theoretical  predictions  using  a  skin  thermal 
model  for  near-IR  laser  pulses  at  1318  nm  and  1540  nm.  For  the  1318-nm  data,  a  Q-switched,  50-ns  pulse 
with  a  spot  size  of  5  mm  was  applied  to  porcine  skin  and  the  damage  thresholds  were  determined  at  1  hour 
and  24  hours  postexposure  using  Probit  analysis  The  same  analysis  was  conducted  for  a  Q-switched,  30- 
ns  pulse  at  1540  nm  with  a  spot  size  of  5  mm.  The  Yucatan  mini-pig  was  used  as  the  skin  model  for  human 
skin  due  to  its  similarity  to  pigmented  human  skin.  The  ED5o  for  these  skin  exposures  at  24  hours 
postexposure  was  10.5  J/cm  for  the  1318-nm  exposures,  and  6.1  J/cm  for  the  1540-nm  exposures.  These 
results  were  compared  to  thermal  model  predictions.  We  show'  that  the  thermal  model  fails  to  account  for 
the  ED5o  values  observed.  A  brief  discussion  of  the  possible  causes  of  this  discrepancy  is  presented,  These 
thresholds  are  also  compared  with  previously  published  skin  minimum  visible  lesion  (MVL)  thresholds  and 
with  the  ANSI  Standard’s  MPE  for  1318-nm  lasers  at  50  ns  and  1540-nm  lasers  at  30  ns. 
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I.  INTRODUCTION 

Recent  developments  in  solid  state  and  chemical  lasers  technology  have  facilitated  the  growing  use  of 
pulsed,  high-peak-power  lasers  at  the  1300-nm  and  1500-nm  wavelengths.  Numerous  medical,  military 
and  industrial  applications  have  been  identified.  While  several  ocular  injury  studies  for  1300  and  1500  nm 
exist,  few  skin  injury  studies  have  been  conducted  at  these  wavelengths,  especially  for  Q-switched 
operation.  Prior  to  the  recent  advancements  in  the  capability  of  these  laser  systems,  particularly  the  in  the 
1540-nm  regime,  little  concern  was  given  to  potential  skin  hazards  associated  with  these  systems. 
However,  higher-energy  systems  at  these  wavelengths  are  becoming  available  in  smaller,  more  reliable 
packages.  Table-top  laser  systems  are  capable  of  producing  microsecond  pulses  with  tens  of  joules  of 
energy  per  pulse  \ 

This  study  provides  the  estimated  dose  (ED50)  required  for  a  50%  probability  of  producing  laser-induced 
lesions  to  skin.  ED^o  data  for  skin  injury  at  1  hour  and  24  hours  postexposure  were  determined  for  high- 
energy,  Q-switched  laser  exposures  at  1318  nm  with  pulse  duration  of  50  ns  and  beam  diameter  of  5  mm. 
In  addition,  ED50  for  1540-nm  exposures  with  pulse  duration  of  30  ns  and  beam  diameter  of  5  mm  were 
determined  at  1  hour  and  24  hour  postexposure.  The  data  were  collected  using  the  Yucatan  mini-pig  ( Sus 
scrofa  domestica)  as  the  in  vivo  model.  The  Yucatan  mini-pig  was  selected  due  to  the  similarity  of  its  flank 
skin  to  that  of  the  human  forearm  skin  The  reported  ED5o  for  each  wavelength  was  compared  to  predicted 
damage  using  thermal  modeling. 
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II.  MATERIALS  AND  METHODS 


For  both  wavelengths  measured,  the  procedure  for  obtaining  the  ED50  was  essentially  the  same.  Laser 
exposures  at  various  energies  were  delivered  to  the  subjects  and  Probit  analysis  was  used  to  determine  the 
ED50.  The  exposures  were  placed  in  the  spaces  of  two  6-cm  x  6-cm  square  grids  on  each  flank.  The  energy 
was  systematically  varied  over  a  wide  range  between  maximum  and  minimum  energy  delivered.  The 
placement  of  the  energy  was  in  a  pseudo-random  order  determined  prior  to  exposure.  The  only  major 
differences  in  the  procedure  for  measuring  the  ED50s  are  in  the  laser  used  and  the  slight  differences  in  the 
optical  setup  for  each  laser  due  to  the  different  methods  for  varying  the  exposure  energy. 

2.1  Animal  prep 

Yucatan  mini-pigs,  weighing  between  15  and  20  kg,  were  used  as  the  in  vivo  model  in  this  study.  Three 
separate  flanks  were  used  at  each  wavelength.  The  study  progressed  under  the  animal  use  protocol  titled 
“Evaluation  of  Laser  induced  Corneal  Lesions  in  the  Dutch  Belted  Rabbit  and  Skm  Lesions  in  the  Yucatan 
Mini-Pig,”  which  was  approved  by  the  Brooks  City -Base,  TX  Institutional  Animal  Care  and  Use 
Committee  (IACUC).  None  of  the  animals  were  euthanized  after  exposure  or  biopsy.  The  subjects  were 
fed  standard  diets  and  had  unlimited  access  to  water.  However,  solid  food  was  withheld  for  12  hours  prior 
to  laser  exposure  and  biopsy  collection.  The  animals  involved  in  this  study  were  procured,  maintained,  and 
used  in  accordance  with  the  Federal  Animal  Welfare  Act  and  the  “Guide  for  the  Care  and  Use  of 
Laboratory  Animals”  prepared  by  the  Institute  of  Laboratory  Animal  Resources  —  National  Research 
Council.  Brooks  City-Base,  TX  has  been  fully  accredited  by  the  Association  for  Assessment  and 
Accreditation  of  Laboratory  Animal  Care,  International  (AAALAC)  since  1967. 

The  Yucatan  mini-pig  subjects  were  received  from  an  attending  veterinary  technician  via  stretcher,  sedated 
with  an  indwelling  intravenous  catheter  placed  in  an  ear  vein,  and  intubated  prior  to  arrival  at  the 
laboratory.  The  pigs  were  sedated  by  a  single  syringe  injection  of  Tiletamine/Zolazepam  (4-6  mg/kg) 
intramuscular  (IM)  and  Xylazine  (2.2  mg/kg)  IM,  and  maintained  on  inhalation  isoflurane  anesthesia 
during  all  procedures.  After  sedation,  hair  on  the  flank  was  clipped  using  hand  clippers  and  the  underlying 
skin  was  cleansed  with  chlorhexadine  solution  and  allowed  to  air  dry.  Prior  to  the  experimental  procedures, 
the  cleansed  skin  was  inspected  by  each  of  three  evaluators  to  check  for  redness,  irritation  or  other 
confounding  marks.  All  subjects  were  kept  warm  dunng  the  entirety  of  the  procedures.  Physiological 
parameters  were  monitored  throughout  all  procedures.  Buprenorphine  (0.05-0.1  mgkg)  was  administered 
intramuscularly  for  analgesia  after  biopsies  were  complete.  The  animals  were  returned  to  their  runs  upon 
recovery  to  sternal  recumbency  from  anesthesia. 

2.2  1318-nm  exposures 

The  exposures  for  the  determination  of  ED50  at  1318  nm  were  delivered  using  a  Nd:YAG,  Q-switched 
laser.  The  laser  produced  a  single  pulse  with  50-ns  duration.  The  energy  of  each  pulse  was  varied  by  a 
polarizing,  beam-splitting  cube  and  X/2  waveplate  combination.  The  cube  and  waveplate  combination  were 
set  up  inside  the  laser  cavity  itself.  The  delivered  energy  of  a  single  pulse  ranged  from  —100  mJ  to  3  J.  The 
delivered  energy'  was  measured  by  placing  a  beamsplitter  in  the  beam  path  and  directing  a  fraction  of  the 
energy  onto  an  energy  probe  (Molectron  JD25  and  JD50),  and  read  with  an  energy  meter  (Molectron 
EPM2000).  This  detector  measured  what  is  known  as  the  reference  energy.  Before  and  after  each 
exposure  set,  a  second  energy  probe  (Ophir  30[150]-A-HE)  was  placed  at  the  location  of  the  subject;  this 
detector  measured  the  delivered  energy.  The  ratio  between  the  reference  energy  and  the  delivered  energy 
was  recorded  for  each  setting  of  the  X/2  waveplate  used  in  the  exposures  (see  Figure  1). 

The  desired  exposure  beam  diameter  was  achieved  by  placing  a  telescope  consisting  of  a  negative  lens 
followed  by  a  positive  lens.  The  spot  size  of  5  mm  was  produced  by  adjusting  the  placement  of  the  subject 
with  respect  to  the  final  positive  lens.  A  co-aligned  argon  laser  was  used  as  a  visible  alignment  beam  to 
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aim  the  IR  beam  onto  the  subject.  A  steering  mirror  was  placed  after  the  telescope  to  make  slight 
adjustments  to  the  vertical  and  horizontal  placement  of  the  alignment  beam.  The  distance  between  the 
subject  and  the  final  positive  lens  was  kept  constant  to  ensure  a  consistent  beam  diameter  for  each 
exposure 
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Figure  1.  1318-nm  experimental  setup.  Exposure  energy  was  varied  by  rotation  of  a  Vfc-wa ve  plate  located  w  ithin  the 
1318-nm  Q-switched  laser. 


2.3  1540-nm  exposures 

A  flashlamp -pumped,  Er:Glass  laser  was  used  to  provide  the  exposures  at  1540  nm.  The  laser  produced 
pulses  in  the  30-40-ns  range  at  various  pulse  energies  up  to  3.5  joules  per  pulse,  The  layout  of  the  optical 
system  consisted  of  a  90/1 0  beamsplitter  followed  by  a  1-m  focal  length  lens.  The  lens  was  used  to  achieve 
the  desired  spot  size  of  5  mm.  A  final  turning  mirror  was  used  to  aim  the  beam  to  the  desired  location. 
Energy  measurements  of  the  laser  light  passing  through  the  beamsplitter  were  made  using  an  Ophir 
Laserstar  energy  meter  with  Ophir  (model  number:  30[1 50J-A-HE)  energy  probes.  A  HeNe  laser  was  used 
as  a  visible  alignment  beam  to  help  locate  the  exposure  point  (see  Figure  2). 

As  with  the  1318-nm  exposures,  energy  measurements  w'ere  made  at  the  subject  location  before  and  after 
each  set  of  exposures.  The  ratio  of  these  two  energy  readings  was  used  to  determine  the  delivered  energy 
for  each  exposure  based  on  the  recorded  energy  through  the  beamsplitter  Due  to  the  limited  spatial  range 
of  the  final  turning  mirror,  the  position  of  the  subject  needed  adjusting  for  some  of  the  exposure  locations. 
A  sling  was  utilized  to  make  the  small  spatial  adjustments.  The  subject’s  elevation  could  be  changed  as 
well  as  the  lateral  position  of  the  subject.  Again,  the  distance  of  the  subject  from  the  laser  was  kept 
constant  to  ensure  uniform  spot  size  for  each  exposure. 
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Figure  2:  1 540-nm  experimental  setup.  Distance  from  the  1 -meter  focal  length  lens  to  the  subject  was  kept  constant. 

2.4  Reading  and  analysis 

During  the  exposures,  notes  were  made  if  an  exposure  caused  an  audible  “p°P”>  a  visible  flash,  or  both. 
The  exposures  were  examined  at  1  hour  and  24  hours  postexposure  Three  trained  readers  independently 
examined  each  exposure  site  and  recorded  either  “lesion’5  or  “no  lesion”  for  each  site.  Once  the  three 
readers  were  finished,  the  individual  results  were  compared.  Agreement  of  at  least  two  of  the  three  readers 
was  needed  for  a  positive  lesion  result  to  be  confirmed. 

The  Probit2  statistical  analysis  package  was  used  to  calculate  the  ED5o  for  each  wavelength.  Supporting  the 
analyses  of  the  ED5o,  the  fiducial  limit  at  the  95%  confidence  level,  along  with  the  slope  and  probability  of 
the  Probit  curve  w^as  also  calculated  at  each  wavelength  at  both  1  hour  and  24  hours  postexposure. 
Following  the  experimental  exposures,  two  5-mm  punch  biopsies  were  obtained  from  each  subject 
immediately  after  the  1-hour  reading  and  two  biopsies  again  after  the  24-hour  reading.  Biopsies  were 
submitted  for  complete  histological  analysis.  A  control  biopsy  w'as  collected  on  a  portion  of  the  flank 
superior  to  any  of  the  exposure  sites.  All  biopsy  sites  were  closed  with  non-absorbable  2-0  sutures  and 
topically  medicated  with  Trio-mycin  ointment  for  infection  prophylaxis.  The  results  of  histological 
analysis  were  not  available  at  the  time  of  publication  and  will  be  reported  in  future  manuscript 


III.  RESULTS 

The  calculated  ED50  for  1318  run  is  based  on  216  exposures  on  three  subjects.  On  each  subject  72 
exposures  were  placed  on  one  flank.  The  1 540-nm  ED50  was  calculated  using  216  exposures  as  well; 
however,  two  subjects  were  used.  Both  flanks  of  one  subject  received  72  exposures  and  one  flank  of  the 
other  subject  was  used.  Table  1  lists  the  results  for  the  1318-nm  and  1 540-nm  exposures  at  both  1  hour  and 
24  hours  postexposure.  Along  with  the  ED50,  the  slope  for  each  Probit  curve  is  given.  The  ED50s  reported 
here  are  above  the  MPE  (maximum  permissible  exposure)  level  as  set  by  ANSI  (Z136  1-2000)  for  skin 
exposures  at  the  reported  wavelengths  and  pulse  durations.  The  MPE  for  the  1318-nm  exposure  condition 
studied  here  is  a  value  of  0.1  Jem'2  The  MPE  for  the  1 540-nm  exposures  is  1.0  Jem'2.  The  results  for  the 
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1 540-nm  exposures  agree  reasonably  well  with  other  work4**  considering  the  difference  in  spot  sizes  and 
pulse  durations. 


The  majority  of  the  lesions  occurred  on  exposure  sites  in  which  an  audible  “pop”  and  visible  flash  were 
observed.  This  suggests  that  the  high  peak-irradiances  incident  on  the  skin  induced  the  formation  of  a 
plasma  at  the  interface  between  the  transport  medium  (air)  and  the  subject  medium  (skin).  This  plasma  was 
facilitated  by  the  presence  of  impurities  (i.e.  dead  skin  cells)  on  the  surface  of  the  skin.  We  will  examine 
the  role  of  this  plasma  as  part  of  the  damage  mechanism  in  the  following  section. 


Experimental  Setup 

Number  of  Subjects  & 
Exposures 

ED50  (J  cm'2) 
1-Hour  Reading 

ED50  (J  cm'2) 
24-Hour  Reading 

Probit  Curve 
Slope  = 
dprob/  Sdose* 

1318  nm 

5-mm  beam  diameter  50  ns 
Three  subjects,  three  flanks 

2 1 6  exposures 

6.5  (7.0 -5.9) 

10.5  (11.1  -  10.0) 

10.7 

1 540  nm 

5-mm  beam  diameter  30  ns 
Two  subjects,  three  flanks 

2 1 6  exposures 

6.3  (6.8-  5.8) 

6.1  (6.5 -5.5) 

6.6 

*  Probit  curve  slope  for  24-hr  ED5(). 

Table  1:  1-hr  and  24-hr  ED50.  Fiducial  limits  are  shown  in  parentheses. 


Figure  4:  13 18-nm  laser  exposures  after  24  hours  showing  slight  reddening  at  lesion  site 
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Figure  6:  1540-nm  laser  exposures  after  24  hours,  showing  lesion  site  compared  to  no-lesion  site. 


IV.  DISCUSSION 


It  is  interesting  to  note  the  differences  in  the  1-hr  and  24-hr  ED50s  for  the  two  wavelengths  .  For  the  1318- 
nm  exposures,  the  ED5q  threshold  increases  from  the  1-hr  to  the  24-hr  exposures.  This  may  indicate  the 
type  of  damage  incurred  by  the  tissue.  It  appeared  that  most  of  the  1  -hr  lesions  were  on  the  surface  of  the 
skin.  The  lesions  were  generally  white  in  color  and  did  not  seem  to  raise  the  skin  in  a  significant  manner. 
Many  of  these  lesions  were  no  longer  visible  at  24  hrs,  suggesting  that  only  the  outermost  layer  of  the  skin 
was  affected.  This  could  possibly  be  due  to  the  pulse  duration  (50  ns)  and  high  peak-irradiances  used  for 
these  exposures.  However,  the  thresholds  for  the  1-hr  and  24-hr  ED50S  for  the  1540-nm  exposures 
remained  virtually  constant.  The  spot  sizes  used  were  comparable,  as  were  the  pulse  durations  (30ns)  and 
peak  irradiances.  This  may  suggest  a  wavelength-dependent  effect  as  to  why  the  lesions  produced  by 
exposure  to  1540-nm  laser  light  possessed  a  higher  probability  of  being  visible  at  24  hrs  post  exposure. 
The  lesions  observed  at  1  hr  for  the  1540-nm  exposures  were  much  redder  in  color  and  raised  the  skin  in  a 
blistering  fashion.  In  addition,  some  of  the  lesions  produced  with  the  1540-nm  laser,  particularly  those 
produced  using  the  maximum  radiant  exposure,  persisted  for  several  weeks,  as  seen  in  the  small  brown 
spots  in  Figure  7. 
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Figure  7:  1540-nm  laser  exposures  after  2  weeks.  The  dark  scar  with  sutures  is  a  biopsy  punch  site. 

4.1  Thermal  damage  modeling 

At  the  wavelengths  studied,  the  normal  damage  mechanism  for  long  pulses  and  continuous  wave  exposures 
is  photothermal.  However,  at  the  pulse  durations  and  energy  levels  reported  in  this  study,  the  exposures  are 
within  the  thermal  confinement  regime9,  suggesting  that  thermal  damage  of  the  tissue  is  unlikely.  To 
accomplish  this  task,  thermal  damage  models  were  used  to  compare  computational  predictions  at  infrared 
w  avelengths  to  physical  data  The  thermal  damage  model  used  in  this  w'ork  is  a  modified  version  of  the 
work  done  by  Henriques1''  The  damage  model  takes  the  basic  form  of  a  rate  equation  for  heat  transfer: 

t 

ri(r,z)=  A  Jexp(“  E  /  RT)dt  (1) 

o 

w'here 

A  =  pre-exponential  factor  (s"‘) 

E  —  activation  energy  for  the  reaction  (caEM) 

R-  2.0  =  universal  gas  constant  (cal/M-K) 

T=  T(r,z,t)  temperature  (K) 

Q(r,z)  -  damage  integral 

t  =  time  at  final  recovery  of  temperature  after  exposure. 

Henriques  set  up  the  rate  equation  such  that  complete  necrosis  of  affected  skin  cells  occurs  when  Q  =  1. 
From  this  condition  the  values  for  A  and  E  are  found  to  be 

A=  3  1  X  1098  (1/s) 

£=  150,000  (cal/M). 

The  Takata  Model  1  used  in  this  study  solves  Eq.  (1)  by  employing  a  “finite  difference”  method  used  by 
Mainster,  ei  al  2  3  It  also  introduces  modifications  to  take  into  account  specific  skin  geometries  as  well  as 
optical  and  thermal  properties  of  the  skin.  Additionally,  the  Takata  model  allows  for  appropriate  alteration 
of  the  beam  profile  as  the  beam  passes  through  various  layers  of  skin. 

Calculations  for  each  wavelength  were  conducted  to  estimate  the  temperature  rise  caused  by  the  laser 
exposure.  The  calculated  temperature  rises  were  then  compared  to  the  minimum  visible  lesion  (MVL) 
threshold  associated  with  temperature  rises  of  skin.  The  results  suggest  that  the  mechanism  for  the  damage 
observed  from  these  laser  exposures  is  not  exclusively  thermal  in  nature.  The  predicted  temperature  rise 
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for  the  laser  exposures  at  1318  nm  were  4°  C  without  modeling  intercellular  scattering  and  8°  C  when  the 
scattering  was  included.  This  is  quite  low  compared  to  the  2nd-degree  bum  MVL  threshold  of  ~29°C.  The 
results  were  similar  for  the  1540-nm  exposures.  A  temperature  rise  of  4  7°C  was  calculated  through  the 
model,  which  again  is  lower  than  the  thermal  ED50  of 29°C. 

4.2  Plasma  formation 

The  model  used  obviously  does  not  account  for  the  damage  produced  by  these  laser  exposures.  While  the 
pulse  duration  and  energy  levels  observed  at  the  reported  ED50  levels  are  withm  the  stress  confinement 
criteria1' ,  the  damage  mechanism  responsible  for  the  lesions  is  not  thought  to  be  ablative  in  nature  The 
ED50  levels  reported  here  appear  to  be  well  below  the  extrapolated  thresholds  determined  from  results  of 
nanosecond  pulses  at  1064  nm.14  As  mentioned  above,  no  histological  results  are  available  at  time  of 
printing  and  these  results  will  allow  a  definitive  answer  on  whether  or  not  ablation  was  occurring  due  to  the 
exposures. 

Although  the  ED50  levels  reported  here  seem  to  be  below  the  thresholds  for  tissue  ablation,  the  levels  are 
capable  of  inducing  a  laser  induced  breakdown  (LIB)  event  at  the  skin’s  surface.  The  ED5o  for  both 
wavelengths  correspond  to  roughly  2x1 08  W/cnV  (2.1xl08  W/cirf  for  1318-nm  pulses,  and  2.03xl08  W/cm2 
for  1540-nm  pulses).  These  irradiance  levels  lie  within  the  plasma  formation  threshold  for  short  laser 
pulses  described  in  Welch9.  Using  an  index  of  refraction^  for  the  skin  of  n~1.4,  the  electric  field  due  the 
laser  pulse  impinging  on  the  skin  for  both  wavelengths  is  -3.3x10  V/m.16  This  result  confirms  the  flashes 
of  white  light  noted  during  exposures  at  and  above  the  ED50  level  were  due  to  LIB.  This  also  explains  the 
acoustic  “pop”  heard  during  the  exposures. 

As  mention  above,  the  1318-nm  and  1540-nm  exposures  are  in  the  stress  confinement  regime.  This 
suggests  that  the  shock  wave  created  by  the  acoustic  pop  is  capable  of  causing  mechanical  damage  to  the 
skin.  Additionally,  the  plasma  itself  may  cause  surface  temperatures  of  the  skin  to  rise  much  higher  than 
the  thermal  model  used  predicts.  Plasma  clouds  of  this  nature  can  reach  temperatures  of  several  thousand 
degrees  K.lr  The  short  time-scales  and  high  peak-irradiances  involved  with  these  laser  exposures  are  not 
characterized  well  with  existing  thermal  models  employed  in  this  work. 


V.  CONCLUSIONS 

ED50  data  are  reported  for  the  wavelengths  of  1318  nm  and  1540  nm  with  pulse  durations  of  50  ns  and  30 
ns  respectively.  Beam  diameters  of  5  mm  where  used  for  both  wavelengths  studied.  The  data  reported  are 
compared  to  ANSI  Z1 36. 1-2000  and  found  to  be  above  the  listed  MPE  for  each  case.  Thermal  modeling 
analysis  is  also  conducted  and  a  discrepancy  is  found  between  the  model  predictions  of  temperature  rise 
and  observed  ED50s.  The  predicted  temperature  rise  from  the  model  is  much  lower  than  the  2nd-degree  bum 
threshold.  Perhaps  the  thermal  rise  of  the  skin  is  greater  than  the  thermal  model  predicts  due  to  the  short 
time  scale  involved.  Moreover,  the  plasma  created  during  the  exposure  may  have  caused  heating  of  the 
tissue  that  the  model  does  not  take  into  account.  Another  possibility  is  that  the  acoustic  wave  created  by 
the  formation  of  the  plasma  caused  the  damage.  More  work  in  the  modeling  of  interactions  between  laser 
radiation  and  in  vivo  subjects  needs  to  be  conducted  for  situations  with  similar  time  scales  and  peak 
irradiances. 
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